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Abstract 



We study the effects of QCD corrections to the process e^e~ — > + X — 
M^uM'u + X above threshold. We show how to treat consistently to 0{as) the 
gluon radiation in both the production and the decay of the top quarks, while 
maintaining all angular correlations in the event. At this order there is an ambiguity 
in the event reconstruction whenever a real gluon occurs in the final state. We study 
the effects of this ambiguity on the top mass and helicity angle distributions. For a 
top mass of 175 GeV and collider energy of 400 GeV the gluon radiation is emitted 
predominantly in the decay of the top quarks. 
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1 Introduction 



Recently, the CDF|T| and the DO|0] collaborations at Fermilab announced the ob- 
servation of the top quark in pp collisions at the Tevatron. Both groups saw a 
statistically significant excess of dilepton and lepton+jets events with the proper 
kinematic properties and bottom quark tags needed to indicate ti production. Fur- 
thermore, they were able to extract mass values for the top quark by fitting to 
events consisting of a single lepton plus four jets. The DO group found a mass 
of 199^^? ± 22 GeV, while CDF obtained a mass of 176 ± 8 ± 10 GeV. Both of 
these mass measurements are in excellent agreement with the value of 175 ± ll^ig 
GeV obtained indirectly from a global fit0 to the electroweak data from LEP and 
SLAG. The direct observation of the top quark at the Tevatron heralds the start 
of a new era in the study of fiavor physics. 

The top quark is certainly unique among the six known quarks. It is by far the 
heaviest; more than 30 times as massive as the bottom quark and even more massive 
than the W and Z bosons. Gorrespondingly, the top quark also has the largest 
coupling to the symmetry breaking sector of all the known particles. This large 
coupling to the Higgs sector may give rise to deviations from its expected behavior, 
thereby offering clues to symmetry breaking, fermion mass generation, quark family 
replication, and other deficiencies of the Standard Model. For example, in top-color 
and extended technicolor (ETG) models the top quark may have non-standard 
couplings to the weak vector bosons 0] or there may be resonant enhancement of 
ti production^. It is of utmost importance to examine the top quark properties 
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as precisely as possible. 

A more basic consequence of the large top quark mass is its short lifetime. For 
large mass the lifetime of the top quark scales as [1.7 GeV ■ (mt/175 GeV)'^]^^, 
and so the top quark decays very rapidly to a bottom quark and a W. Thus, 
unlike the lighter quarks which form hadronic bound states before decaying, the 
top quark behaves more like a heavy lepton, decaying as an unbound fermion. In 
fact, it decays long before depolarization, so that its spin information can be easily 
reconstructed from the momenta of its decay products. This fact will be extremely 
useful for extracting information about the top quark parameters. 

An ideal place to study the top quark is in e'^e~ collisions!^, |^, where the 
colorless initial state provides a clean event environment, and there is the possibility 
of initial-state polarization. By varying the beam energy it is possible to scan the 
threshold region or to study the top above threshold. There have been many studies 
of top production near threshold, where the resonance behavior can be calculated 
in perturbative QCD and the top mass can be obtained to high accuracy ||]. In 
this paper we will instead concentrate on the continuum tt production. At tree 
level the event is characterized by six final-state particles arising from the process 
e~^e^ ^ tt ^ bW~^bW~ —>■ M'^ubi"!?. These six particles contain a wealth of 
information in their relative momenta, angles, and polarizations. By reconstructing 
the helicity angles of the top quarks and the W^s, it is straightforward to extract 
the top quark parameters. 

Although the top quark is produced and decays essentially as an unbound 



2 



fermion, it still feels the strong interactions and will radiate gluons — both in its 
production phase and its decay phase. Thus, it is useful to see how the tree- 
level picture and experimental analysis will be affected by QCD corrections. The 
0{as) corrections to the production have been studied in several papers, including 
analyses of the effects on production angle distributions ^ and polarizations |[T0[| . 
Similarly, studies of the 0{as) corrections on the top decay have been done, with 



analyses of energy distributions, and angular distributions from polarized tops|[TT 
However, the top production and decays do not occur in isolation from each other. 
For events with an extra gluon jet it is not a priori obvious whether to assign the 
extra jet to the production, the t decay, or the i decay. At the very least, the extra 
jets will add one more degree of complexity to the event reconstruction process. 
Therefore, it is necessary to assess the impact of these radiative corrections on the 



full event f 12 



To this end we have constructed a next-to-leading order (NLO) Monte Carlo 
which treats consistently to 0{as) the radiative corrections to both production 
and decay of the top quarks. To set the stage for this NLO analysis we begin 
by reviewing the salient features of the e~^e^ tt event at tree level using helic- 
ity decomposition in section 2. Then in section 3 we analyze the cross section at 
next-to-leading order and give the details of the Monte Carlo, describing the ap- 
proximations used and the methods for subtracting the infrared (IR) divergences in 
production and decay. We also include two appendices with the helicity amplitudes 
for top production and decay with real gluons. In section 4 we use the Monte Carlo 
to study the effects of gluon radiation on the top quark mass measurement and to 
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re-examine the helicity angle distributions at next-to-leading order. In this section 
we assume that the only ambiguities are in the placement of the extra gluon jet, 
that the VF's and bottom quarks are correctly identified, and we investigate how 
the distributions vary with the algorithm used for assigning the gluon jet. Then in 
section 5 we make another pass through the mass distributions with more realistic 
experimental assumptions for the event. The purpose of this section is to identify 
which physical inputs have the largest effect on the continuum measurement of the 
top mass. In section 6 we offer our conclusions. 

2 Review of the tree-level analysis 

Even at tree level the full e^e~ tt event is quite complex. The six-particle 
final state can be characterized in many possible ways by the relative momenta 
and angles in the event. It is an important conceptual problem to clarify which 
pieces of information are most important, and how all of the various kinematic 
measurements available cooperate to illuminate the basic physics. The solution to 
this problem is suggested by the fact that the event is actually a series of on-shell 
two-body decays: j*,Z* tt, t ^ bW'^, and i'^u. Thus, by considering 

intermediate states of definite helicities, the event is highly constrained simply by 
conservation of angular momentum. The different helicity states are revealed by 
the angular distributions of their decay products, while the relative amplitudes for 
the different helicity combinations are easily related to the couplings at the top 
quark production and decay vertices. In this section we describe this tree-level 
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helicity analysis. Although this has been discussed before in the literature, most 



notably by Kane, Ladinsky, and Yuan|13|, we will review it here for pedagogical 



purposes and to set the notation for the discussion of QCD corrections. 

The dominant effects of new physics on the process e^e~ ^ tt ^ bW~^bW~ 
can be described in terms of form factors included at the production and decay 
vertices. The t — > bW~^ decay vertex can be written 

^M^^ = ^j=^{nFr,PL + + ^[^^T^^ + F^^rPl]} , (1) 

where = (1 ± 75)/2, and we have neglected a third pair of form factors which 
do not contribute to decays to on-shell W^s or massless fermions. We have chosen 
the subscripts L, R of the form factors so that they indicate the helicity of the 
outgoing bottom quark in the limit mf, = 0, which we will use in all of our matrix 
element calculations. At tree level in the standard model = 1 and all other 
form factors are zero. In fact, = F^ = to all orders in the standard model 
in the limit of massless bottom quark. The antitop form factors are identical to 
these in the limit of CP invariance. 

Similarly, the 'j, Z ti production vertices can be written 

^M^^ = te{r[Fly + Fl^^,] + ^^[F^v + F^^^,]} , (2) 

where each form factor can be a function of the center-of-mass energy ^/s, the 
superscript is i = 7, Z, and we have again dropped a third pair of form factors 
which are unobservable. At tree level in the standard model F^y = |, Ffy = 
(1 ~ '^sIj)/ SwCw, and Ff^ = {—j)/swCw, and all others are zero. Here, = sm9w 



and Cyj = cosOw In the limit of CP invariance = 0. The production analysis 
is simplified if we consider separately the two possible helicities of the incoming 
electrons, so that the contribution of the photon and the Z add coherently. We 
define new form factors by 



where the subscripts, i — 1,2 and j — V,A refer to the structure of the form factor, 
and the superscripts refer to the helicity of the incoming electron. 

We are now ready to discuss the helicity angle description of the complete event. 
As mentioned previously, in the limit of narrow width for the top and the W, the 
event can be considered as a succession of two-body decays. The first process we 
consider is the decay of the virtual 7, Z boson into the tt pair. Note that the 
intermediate vector boson receives twice the helicity of the initial electron, along 
the beam direction. This process can be described in the e'^e~ center-of-momentum 
frame by two angles, the polar angle 9 and the azimuthal angle 4> of the top with 
respect to the electron beam axis. Using the notation and tn to denote the 
helicities ht — —1/2 and ht — +1/2, we obtain the matrix elements 




(3) 



M{eLeR 



M{eLeR 



M{eLeR 



tRtR) 



tRh) 



tLtR) 



tLh) 



[^w - P^iA + ^2v\ (1 + COS d)e-"^ 
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MieRCL^tLtR) = -[J^i^^-/3J^5i + J^2y](l-cos^)e^'^ (4) 

MieRCL^tRh) = -[J^^v + P^^A + ^2v]i'i- + cose)e"f' 

M^eReL^tLh) = r\^w + l\^2v + (^^2A)]{^^^0)e"f> 

MicRCL ^ tRtR) = r^[T^v + l\^2v - /9-^2a)] {si^e)e"^ , 

where we have removed a factor of ie^. Here, (3^ — {1 —4m^/s) and 7 = ^/s/ (2mt). 
For longitudinally polarized beams the dependence will vanish. 

The nice aspect of this helicity formalism is that the angular dependence of 
each of the amplitudes is determined, up to a relative phase, simply by angular 
momentum conservation. For instance, in the first matrix element the virtual vector 
boson has helicity -1 along the electron beam direction, the top has helicity -1/2, 
and the antitop has helicity +1/2. To conserve angular momentum the top must 
move in the electron direction and the antitop must move in the positron direction; 
hence the (1 + cos^^) dependence. By measuring the angular distributions it is 
straightforward to extract the relative weights for each helicity combination, and 
thereby obtain the top quark form factors. 

As an example, we plot in Fig. 1 the tree-level Standard Model production cross- 
section as a function of cos 9 for a top mass of 175 GeV and a coUider energy of 400 
GeV for polarized electron beams. We have also plotted the hehcity subprocesses. 
Here we see that the e^'s produce predominantly t^s highly peaked in the forward 
direction, while e^'s produce predominantly Irs peaked in the forward direction. 
This can easily be understood in the limit of high energy, where the SU (2)^ x C/(l) 



7 



symmetry is restored and the squared matrix elements become 
\M{eLeR 
\M{eLeR 
\M{eReL 
\M{eReL 

while the remaining matrix elements vanish. Thus, longitudinally polarized elec- 
trons are an excellent source of polarized top quarks. 

The next stage in the event is the decay of the top t bW^. This process 
is most conveniently described in the top rest frame obtained from the lab frame 
by rotating the axes —(f), then —9, and then boosting in the direction opposite to 
the top momentum. The helicity angles in this frame are the polar angle Xt a-nd 
the azimuthal angle ipt of the W boson with respect to the top momentum axis. 
Using the notation (L, R, Z) to denote the helicities ( — 1, +1, 0), we obtain the 
helicity amplitudes for the left-handed bottom quarks: 

MitR^h^W^) = w-'[Fr,-\w'Fl]{cos^)e^^^'' 

MitR^bLW^) = V2[F,^-lF2^](-sin|)e^'^*/2 (6) 
Mit,^b,W^) = v^[F,^ - If^^] (cos |)e-^*/2 
M{tL^bLW+) = M{tR^bLW+) = 0, 

where w = mw/mt, and we have dropped an overall factor of igmt{l — w'^y^'^/\/2. 
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Mr) I' 



(-L + ^)^(l + cos^f - 1.41(l + cos^)^ 



1 \2 



tRtL)^ = (tTt) (1-COS^)' ~ 0.19 (1- COS 



tRtL)\'' = (-|-)'(l + COSef ~ 0.75 (1 + COS 



(5) 



tLtR)\' 



cost 



0.05(1 - COS , 



The matrix elements for right-handed bottom quarks are obtained from these by 
replacing everywhere L R, ipt — and Xt ^ ~Xt- 

As before, the angular dependence is exactly what is expected from angular 
momentum conservation in the decay of a spin- 1/2 object. In addition, in the 
Standard Model in the limit irib — the top can only decay to 6l's. Therefore, 
it must decay to Wz"'s in the direction of the top quark spin, to W^/'s in the 
direction opposite to the top quark spin, and it cannot decay to W^^s at all. In 
Fig. 2 we display this by plotting the tu — > bW'^ decay distribution as a function 
of cosxt, while also plotting the helicity subprocesses. For increasing top mass the 
distribution becomes more sloped in the forward direction, indicating an increased 
partial width to W^- 

The antitop decay i hW^ can be described in an analogous manner in the 
antitop rest-frame, obtained from the lab frame by rotating the axes —0, then 
TT — 9, and then boosting in the direction opposite to the antitop momentum. The 
helicity angles in this frame are the polar angles of the W~, Xt and 'ipt, with respect 
to the antitop momentum axis. If CP is a good symmetry we can obtain the matrix 
elements using 

M{tn ^ b,W+) = M{t^h ^ KWlx) , (7) 
while replacing xt Xt and ipt —>■ —ipt- 

The final step in the decay chain is W'^ iiy. We work in the W'^ rest frame 
obtained from the top rest frame by rotating the axes —il^t, then — Xt, and then 
boosting against the momentum. The helicity angles in this frame are the 



9 



polar angle x and the azimuthal angle of the charged lepton with respect to 
the W'^ momentum axis. For hadronic decays we can just replace i'^ with the 
antiquark and u with the quark. In the Standard model the W'^ can only decay to 
i'jii'L in the limit of massless leptons. The helicity amplitudes are 

M{W^^i+u) = -L(i + cosx)e^^ 

M{W^^£+u) = sinx (8) 
M{W^^£+u) = -L(l-cosx)e-^^ 

where we have removed a factor of igruw/V^- In Fig. 3 we plot the cosx distribu- 
tion in the W'^ — > i'^u decay, along with hehcity subprocesses, for W'^ produce in 
top decays. The zero at cosx = 1 indicates the absence of right-handed W'^'s. 

Lastly, the decay W~ £~P can be described in the W~ rest-frame, obtained 
from the top rest frame by rotating the axes —ipt, then — Xt, and then boosting 
against the W~ momentum. The helicity angles in this frame are the polar angles 
of the negatively-charged lepton, x and ■0, with respect to the W~ momentum axis. 
We can obtain the helicity amplitudes from 

M{W+ i+u) = M{WZ^ rv) , (9) 

while replacing X ^nd — > —'4). 

In practice, in order to optimize the analysis of the top quark form factors it is 
necessary to study the event in a multi-dimensional space of all these angles. The 
use of helicity angles makes it easy to discern which variables are most important 
for studying which form factors. For example, by cutting on the production angle 
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6, while using a polarized electron beam, it is possible to obtain a sample of highly 
polarized top quarks. With these, one can study the decay form factors by looking 
at both the top decay angle Xt and the W decay angle in order to determine 
the helicities of the W^s. Perhaps the optimum technique would be to use all of 
the helicity angle information in a maximum likelihood fitfl^. In any case we now 



obtain the full tree-level correlation information of the event from 

hh'XX'pp' 

X M{th' bpWy)M{Wy rZ/)p (10) 

for each initial-state helicity configuration. 



3 The event at 

In the narrow top-width approximation, in which the top quarks are treated as on- 
shell particles in the matrix elements, the 0{as) corrections can be unambiguously 
assigned to the ti production process, or to the t-decay or t-decay processes. We 
have constructed a NLO Monte Carlo by separately building a generator for tt 
events with an extra gluon in the production, the t-decay, and the t-decay processes, 
as well as for events with no extra visible gluon. To see how this is implemented it 
is easiest to temporarily ignore the angular correlations. Then the total differential 
cross section, datot, for the event e+e~ ^ ti + X ^ M^uM^u -|- X is just the 
product of the ti + X differential production cross section, da, times the t and i 
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decay distributions: 



To 0{as) this can be written 



w ^ ^^^^ (^^^ 

dcTtot = da . (11) 



d<jlt = -^p^^ + da + da 

, ndV^dV^ 2T^ , ndT^dV^ , , 

+ " T^^^ "(f^' ^^^^ 

where the first term is the tree-level event, the second term includes 0{ol^ correc- 
tions to the tt production, the third and fourth terms contain the corrections to 
the i and the i decay respectively, and the last term is the 0{a.^ correction to the 
widths in the denominator. Note that on integrating over the decay phase space, 
the last three terms cancel so that af^"^^ — a^ + a^ ^ i.e., the integrated total event 
cross section is not affected by the corrections to the top quark decay, as required. 

The 0{as) corrections to the production and decay can be separated into three 
pieces — the virtual (v), soft-gluon (s), and real-gluon (r) contributions: 

da^ = da"" + da'{xo) + da''{xo) 

dr^ = + dr{yo,zo) + dr{yo,Zo) . (13) 

The arbitrary distinction between "soft" and "real" gluons is implemented using 
artifical cutoffs Xq, t/o? -^o? which we will describe more fully below. The real gluons 
are defined to be those produced above the cutoffs and are treated using the exact 
3-body phase space. The soft gluons are those produced below the cutoffs and 
are integrated out analytically, leaving an effective 2-body phase space. Both the 
virtual and the soft contributions are infrared divergent, but their sum is infrared 
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finite. Thus, we can combine the virtual and soft contributions, and we can con- 
veniently separate the full 0{as) cross section into the sum of four sub-event cross 
sections: 



(r 



0^2 



+ Jfof + (f^p • (14) 

The last three contributions have 7 final-state partons, containing a real gluon in 
the production, the t decay, or the i decay respectively. Each of these terms is 
manifestly positive-definite. The first contribution has only 6 final-state partons 
and is given by the following sum: 

datr\x.,y.,z,) = {l-^-^)da'^^;^ + da('"^^\xo)^^;^ (15) 

o rfr(-+-)(yo,;^o)dr° , dTUT^^+^\yo,Zo) 
+ ^2 + ^2 • 

This term may be negative for small values of the cutoffs. A separate Monte 
Carlo is used to generate events for each of the four terms in (0) with all angular 
correlations included. 

We now elaborate on the infrared cancellations, as well as the separation into 



"soft" and "real" gluons, that are used in equation ([I^). The virtual corrections 
to the production and decay processes can be written as corrections to the form 
factors (0) and (Q), with the understanding that they are only expanded to 0{as) 
in the squared amplitudes (PH]). Using dimensional regularization with D = 4 — 2e 
we obtain the correction to the production form factors: 

Zn 
13 



2V 



ZTT 



(16) 



where 
h 



mr 



ra + e) 



yin — — - + ZTT 



+ 



2(3 



2(3 
3 



+ /5 



3 Vi + /3y^2V 1 



it: 



(3^2 



1-/3 



4/3 



In- 



1 + /3 



ZTT 



(17) 



a = 4/3, and /7 



g — t:/'-', "ii^ ~ 3' /-A ~ 0' /y ~ (4 3S^)/s^c^, and — ( 4)/"Su;C^ 



are the tree- level couplings. Also, we have used the Spence function Li2{z) = 
— Jq dtln{l — t)/t. This agrees with the previous results given in Ref. 0. Note 
that the contribution from Re/i is proportional to the tree-level cross section, while 
ImJi does not contribute at 0{as)- 

For the real gluon corrections to ti production it is convenient to define the 
gluon phase space in terms of the variables 



X 



Tmax 



COS0*)/2 



'1^ 



nmax 
'9 



where the maximum energy of the production gluon in the lab frame is 
/3^-\/i/2, and 9^^ is the angle between the gluon and top quark momenta in the ti 
rest frame. The full phase space isO<x<l,0<A<l with the soft gluon 
limit given by x — > 0. Integrating out the gluons in the region x < Xq, for small 



14 



Xo, it is possible to absorb this soft-gluon contribution into the form factors ( [T^ ) 



by replacing Ji ^ Ji + / 



(soft) 

1 ) 



with 



(soft) 



r(l + e) <( ( - -21nxo 



l + l±i?ln^-^ 

2/3 



1 + 



In 



+ 



2/3 



4/34 /3 1 + /3 
ln/32 In 



1 + /32 r 1 - /? 



1 + /? 3 
l-/3x 1/, 1 



(19) 



+2Li,(i^) + ^(lni^^' 



The sum of the virtual and soft contributions Re(/i + /-[''°^*-') is now IR finite. The 
"real" gluons with x > xq are treated using exact kinematics. The matrix elements 
can be written in terms of helicity amplitudes as in section 2. We leave the details 
of this to the appendix A. 

The virtual corrections to the top decay form factors at 0{as) are 



w 

IL 



agCq ( 47r/i 



27r V 



r(i + e) + - 



Infl - 



1 r 5 



^ +ln(l -w;2) 
4 



+ -ln(l -w^) -Li2(w^l 



(20) 



-w m l — W ) 

27r ^ ' 



For the phase space of the real gluon in top decay we use the variables 



max 
9 ' 



cos^,;)/2, 



(21) 



where the maximum energy of the decay gluon in the top rest frame is 



max 

a 



rrit /2){1 — w"^), and is the angle between the gluon and bottom quark momenta 
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in the bW rest frame. The gluon becomes soft in the hmit y and colhnear 
in the hmit z ^ 0. Integrating out the soft and coUinear gluons for which y < yo 
and/or z < Zq, for smaU yo,zo, we can absorb these contributions into the form 
factor F^. They contribute 



271 \ mf J 2e^ e 



--\n(l-w^) 
4 ^ ' 



+4 + ^7^ ^ + 0^1- 

8(1 — w^) ^ 

-^ln(l-w;2) + Li2(l-w;') (22) 
Inw 



4(1 -w;2)^ 
-(1 + lnxo)(l + Inzo) + ^lii-2o| , 

so that the sum of the virtual and soft-gluon contributions SF^ + SF^^^°^^^ is IR 
finite. As in the production process, the "real" gluons with y > yo and z > zo are 
treated using exact kinematics. The helicity amplitudes are given in the appendix 
B. 

It is useful at this stage to describe the Monte Carlo more fully. It is written in 
the C++ programming language and contains a separate event-generator class for 
each of the four sub-channel processes in equation ([I^) . Each of these sub-channel 



generators are in turn derived from a single tree-level generator which produces the 
helicity angles of the event with the exact tree-level distributions. The sub-channel 
generators then produce the relevant gluon kinematic variables, prepare the particle 
four- vectors, and give the event a weight. The production-gluon class generates the 
gluon variables (0) with a soft gluon distribution, while the decay-gluon classes 

16 



generate the gluon variables (^) with a soft and colhnear gluon distribution. This 
results in a very efficient Monte Carlo for each of the four sub-channels. 

The relative contribution from the four sub-channels depends on the artificial 
IR cutoffs (xo, Uo, Zq). The choice of values for these parameters is determined by 
several considerations. First, the analytic integrations of the soft gluons contained 
in (rU) and P^ ) are valid up to terms linear in the cutoffs, so they should be 



kept as small as possible. In addition, they should lie below any physical cutoff, 
determined by the detector energy resolution or the jet definition. However, for very 
small cutoffs the contribution containing the virtual and soft gluons will become 
very large and negative, and there will be large cancellations between it and the 
other sub-channels. Thus, the cutoffs should not be too small or else the numerical 
errors will become prohibitive. Luckily, this last constraint turns out to be not too 
restrictive for our Monte Carlo. For each plot in the next two sections we have 
checked that the results do not change significantly for smaller values of the cutoffs. 
As a final test of our confidence, we have checked that our Monte Carlo reproduces 
the 0{as) production^] and decay distributions [|ll| of previous analyses. 

Our Monte Carlo also allows the inclusion of width effects by generating Breit- 
Wigner resonance distributions for the tops and the W^s. In addition, the kine- 
matic effects of the bottom quark mass can be included. Momentum conservation 
is maintained by shifting the energies of the final-state particles, while keeping the 



helicity angles and the gluon kinematic variables ( |T8| ) and (^if ) fixed. This pro- 
cedure should be good to OlTt/mt) except very near threshold. Note, however, 
that the matrix elements, and hence the event weights, are always computed in the 
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zero-width and nib = limits. Finally, initial state radiation (ISR) can be included 
by generating electron and positron momentum fractions z with the distribution 



function given by Fadin and Kuraev[15 



D,{z) = /3/2(l - (1 + 3/3/8) - + z)/4 , (23) 

where /3 = {2a /n) (Ins /ml — 1). 

It must be noted that the narrow-width approximation is necessary for the 
NLO analysis of this section. As a consequence, the Monte Carlo does not include 
the effects of interference between gluons emitted in the production and gluons 
emitted in the decay. These perturbative effects have been studied in the soft- 
gluon limit in Refs. ||16|. Typically, the interference is only important for gluons 
with energy Eg <Tt. However, it should be considered in a complete analysis. In 
addition, because the final-state bottom quarks do carry bare color, there will be 



some nonperturbative information connecting them in the form of soft hadrons||17||. 
We have neglected this effect here. 



4 The effects of radiated gluons 

In this section we will study the top quark mass reconstruction and helicity angle 
distributions at next-to-leading order. We do this by starting with an ideal event 
situation — no ISR, an ideal 47r detector, perfect partonic-level particle identifica- 
tion. In the subsequent section we will make each of these factors more realistic 
experimentally. The purpose here is to develop our intuition by isolating the purely 
theoretical QCD effects at NLO. If we assume that both bottom quarks and W^s 
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are identified and signed and that there is An detector coverage, then the only 
ambiguity is in where to put the gluon. Does it belong to the t, to the i, or to 
neither? 

Here, we will make this assignment of the real gluon in analogy with the typical 
jet-clustering algorithm used at e+e~ colliders. Defining the quantities /i^ = {pb + 
PgY and /i^ = {pb + PgY, we make the assignment: 

if fi < Ji and /i < ficut =^ gluon belongs to t decay 

(Pt =Pg+Pb+ PW+ , Pt=Pb+ PW-) , 

ii Jx < 11 and /2 < Hcut =^ gluon belongs to i decay 

(24) 

{pt =Pb+ PW+ , Pt=Pg+Pb + Pw-) , 

else ^ gluon belongs to production 

(Pt = Pb + Pw+ , Pt = Pb + Pw-) ■ 
In the limit m;, = 0, we recognize ficut as an infrared cutoff on both the collinear 
and soft gluons in the event. In fact, we can consider the decay gluons to be 
clustered with the bottom quarks|]18[ using the standard jet resolution parameter 



Vcut = f^cut/ ^- varying ficut we change the fraction of events with gluons that are 
not combined with the b or b, and thus are considered to be part of the production 
process. This fraction is plotted in Fig. 4 as a function of the center-of-mass energy 
for various values of ficut- At this fixed order in perturbation theory, the fraction 
can be greater than one, indicating that a resummation of the large logarithms in 
Hcut or fil^^/mf is necessary. As in all of our plots we use a standard top mass of 
175 GeV and a, = 0.12. 

We now consider the top quark mass distribution at ^/s = 400 GeV. Using 
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the algorithm (0), each event produces two mass values = pf and = pf, 
which are binned independently. To see most clearly how the radiation affects this 
distribution, we plot it in Fig. 5(a) in the strict zero top-width limit for values of 
fJ'cut = 5, 10, 20, and oo GeV. Note that for ^cut = oo, all of the observed gluons 
are assigned either to top decay or antitop decay, and none to the production. The 
Monte Carlo cutoffs used are xq = 0.02, yo = 0.005, and zq = 0.01. The 5-function 
spike in the central bin arises from those events in which the top momentum is de- 
termined correctly from its true decay products. The excess below the 5-function 
corresponds to events where a decay gluon is assigned incorrectly and is not in- 
cluded in the top momentum reconstruction. These missed-gluon events become 
less likely as ficut increases, but even for ficut = oo there is a remnant of events 
where the gluon gets assigned to the wrong-charge top quark. The excess above 
the (5-function corresponds to events where an extra gluon is incorrectly included 
in a top momentum reconstruction. This region has two separate contributions, 
from mis-assigned decay gluons and from mis-assigned production gluons. Both 
of these increase with increasing ficut, with the production gluons adding a second 
bump for larger values of this parameter. The deficits in the distribution directly 
on each side of the spike are due to the artificial cutoffs xq, yo, and zo- 

The 5-function peak in this distribution is an artifact of the zero-width ap- 
proximation. Turning on the Breit-Wigner resonance for the top quark effectively 
smears over the 5-function and results in a well-defined IR-finite mass distribution. 
In Fig. 5(b) we plot this distribution using the same values of ficut as before. For 
comparison, we also plot the initial Breit-Wigner distribution. We now choose the 
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Monte Carlo cutoffs to be xq = 0.002, yo = 0.0013, and zq = 0.0028. Tliese cut- 
offs ensure tfiat all production gluons with Eg > 100 MeV and all decay gluons 
with ii,p, > 5 GeV are treated with exact kinematics. The distributions do not 
change significantly for smaller xo,yo,zo- For ficut = 5 GeV we see that the mass 
distribution is severely distorted, while for higher values of Hcut it quickly regains 
an approximate Breit-Wigner shape, with a small decrease in the peak and an in- 
crease in the tail regions. We cannot take the /icut — 5 GeV curve too seriously, 
however, because for small values of ficut "we are probing the collincar-gluon region 
of the decay phase space. On the other hand, the effects of soft-gluon singularities 
are inconsequential, because soft gluons have Eg Ki and do not affect the mass 
measurement. For Hcut ^20 GeV these perturbative mass distributions should be 
reliable. Fig. 5(b) suggests that perhaps the best approach to mass reconstruction 
at = 400 GeV is to treat each extra gluon as coming from decay, combining 
it with whichever top quark has the smaller value of fi. This is because 400 GeV 
is still not too far from threshold, where real gluon radiation in the production 
process is suppressed. 

At higher energies the situation changes dramatically. In Fig. 6 we plot the 
mass distributions at ^/s = 1 TcV for ficut = 5, 20, 80 and oo GeV. At this center- 
of-mass energy we choose Xq — .0001 so that production gluons with Eg > 100 Mev 
are treated with exact kinematics. The best resonant peak occurs for /icut ~ 20 
GeV. At this high energy there is substantial coUinear radiation in the tt production 
process, so that for larger values of Hcut an extra gluon is usually included with 
one of the tops, resulting in a too-large mass reconstruction. These curves are 
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suggestive of the degradation that will occur at this energy, but a resummation of 
the collinear gluons would be necessary to obtain an exact prediction. Certainly, 
determining the top mass at i/i = 1 TeV would be more difficult than at lower 
energies. 

We now turn to the top production angle distribution. For the remainder of this 
section, we work in the strict zero-width and m?, = limits. The production angle 
distribution has been studied before at 0{as) for the pure ti production process 
in 1^. Here we include the effects of radiative corrections in both production and 
decay of the top quarks. Although the corrections to the decay process do not 
affect this distribution for perfectly reconstructed tt events, they are significant 
when reconstruction ambiguities are considered. For a given value of ficut we can 
use the algorithm (^) to reconstruct each event and then bin with respect to the top 
and antitop production variables cos 6 and — cos 6. The tree-level production angle 
distributions for = 175 GeV and ^/s = 400 GeV were shown in Fig. 1. In Fig. 7 
we plot the deviations from the tree-level distribution for several different values of 
/icut for left- and right-handed electron beams. We also plot the pure production 
corrections [^], which assume perfect gluon discrimination and event reconstruction. 
For both electron polarizations the 0{as) corrections tend to increase the slope of 
the distribution with production angle. However, the treatment of the radiative 
gluon can have a significant effect on this correction. For a left-polarized electron 
beam, using smaller values of ficut, the correction even changes sign. This is shown 
further in Table 1, where we give the 0{as) corrections to the forward-backward 
asymmetry of the top quarks for the different values of ficut- 
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In Fig. 8 we examine the effects of tlie gluon ambiguity on the decay angle of 
the top quark to the boson, xt- Using the algorithm (^) the boson is re- 
constructed correctly, but the observed momentum of the top quark, and therefore 
the observed value of Xt, is affected by the treatment of the radiative gluon. In 
Fig. 8 we plot the fraction of observed values of cosxt falling in each 0.1- width bin 
for events with true values of cosxt between -0.1 and 0.0. For small ficut the recon- 
structed values of cos Xt tend to be larger than the true values. The missed gluons 
in the decay lead to an underestimate of the top momentum, which results in an 
underestimate of the angle between the and the top momenta after boosting to 
the top rest frame. As in the previous examples, the most accurate reconstruction 
occurs for large ficut- 

5 More detailed analysis of top mass reconstruc- 
tion 

In this section we re-examine the top mass distribution with more realistic ex- 
perimental assumptions. The neutrinos are undetected and the quark jets are 
indistinguishable. We include the effects of initial-state radiation, and we impose 
simple lab-frame angular cuts to approximate the effects of the detector. We also 
examine the effects of parton energy smearing due to the detector resolution. How- 
ever, we stop short of including final-state hadronization. This analysis is strictly 
at the partonic level. 

We will consider the reconstruction of the top quarks in both the lepton-|-jets 
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mode and the all-jets mode. We require that all of the visible partons must satisfy 



I cos 9iab\ < 0.9, and we cluster |18| the colored partons into jets using the jet reso- 
lution parameter ycut = /^Lt/ ^ with ficut = 20 GeV. We do not consider the effects 
of 6-tagging, treating all hadronic jets as indistinguishable. We then use a simple 
algorithm for tt event reconstruction in each mode. Certainly, these methods can 
be improved and optimized, but they will be sufficient for our purposes. 

In the all-jets mode we require that there be > 6 jets after the cuts and the 
clustering. We then choose two pairs of jets to form the W^s by minimizing the 
quantity 

(Pi + P2)'^ - m'^w + iP3+P4)^-m'^ (25) 

over all combinations of jets. We then combine one or more of the remaining jets 
with each of the W^s, so as to minimize the mass difference between the resulting 
top quarks. 

In the lepton+jets mode we require that there be a charged lepton and > 4 jets 
after the cuts and clustering. The neutrino four-momentum is defined to be equal 
to the missing momentum in the event, p^^ = ptotai — J2Pvisibie, with the additional 
requirement that 

|m(£i/) -mv^l < 10 GeV . (26) 

Then a pair of jets is chosen to form the second W boson by minimizing \{pi + 
P2y —fn'^l over all of the jets. Finally, we combine at least one of the remaining jets 
with each of the W^s, so as to minimize the mass difference between the resulting 
top quarks. 
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We begin our study by including the initial-state radiation, but omitting the 
final-state energy smearing. The mass distributions for the all-jet channel and 
for the lepton+jets channel are shown in Figs. 9(a) and 9(b), respectively, for 
nit = 175 GeV and ^/s = 400 GeV. For comparison we also show the original 
Breit-Wigner distributions, as well as the mass reconstructions at tree-level for both 
channels. The 0{as) distributions exhibit a moderate degradation as compared to 
tree-level and also as compared to the ficut = oo curve of Fig. 5(b) from the previous 
section. This is due to the additional complexity in clustering the radiated gluon 
and reconstructing the event. Naturally, these effects are more serious in the all- 
jet channel. In the lepton+jets channel there can also be errors in the neutrino 
reconstruction due to initial-state radiation. This is the source of the enhanced tail 
at higher masses. Of the all-jet events, 41% survive the cuts and are identified as 
a 6-jet event, while 4% are identified with 7 jets. Of the lepton+jet events, 35% 
survive the cuts and are identified with 4 jets, while 7% are identified with 5 jets. 

In Fig. 10 we show the same distributions with the final-state partons smeared 
in energy to approximate the effects of the detector energy resolution. The hadronic 
and leptonic final-state partons are gaussian-smeared with the parameters used in 
the JLC studyig]: 



had n A J^'P , , , , 

(27) 



af^ 0.4 0.15 



E ' E ' 

where E is in GeV. The smearing has no effect on the efficiency in the all-jets mode, 
but it does reduce the efficiencies in the lepton+jet modes to 22% (4 jets) and 4% 
(5 jets). This is because, when the jet energies are smeared, the reconstructed 
neutrino is less likely to meet the constraint (^^j- From Fig. 10 we conclude that 
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the major contribution to the error on the top mass distribution will probably 
come from the detector energy resolution, making a direct width measurement 
virtually impossible. The gluon radiation also contributes a significant amount to 
the widening of the peak, especially in the all-jets reconstruction channel. As we 
have shown in this paper, this QCD radiative contribution is directly calculable in 
perturbation theory. 

The plots in this section are representative of the accuracy that may be obtain- 
able in a direct mass measurement, although certainly the reconstruction algorithm 
can be better optimized, and 6-tagging would be very useful in this regard. As for 
the angular distributions, we would expect the detector resohition effects to be less 
serious because detector angular resolution is usually better than energy resolution. 
However, the reconstruction errors may still be significant for these distributions. 



6 Conclusions 



As in any hard scattering process, the e"^e~ — > tt event is certainly more complex 
than the basic tree-level parton cross section would indicate. The first step to a 
more realistic treatment should include QCD radiation in the final state. This 
requires the correct handling of radiation both in the 7* tt production process 
and in the t — > bW'^ decay process. In this paper we have shown how to include this 
radiation to 0{as) and have constructed a Monte Carlo generator to study these 
effects. In doing this we have made strong use of the helicity angle formalism, 
which is the most natural for investigating the properties of the top quark. 
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The treatment of the ti event at 0{as) introduces reconstruction ambiguities 
whenever there is real gluon radiation. We have shown how this can alter the top 
mass distribution and the angular distributions. By including the Breit-Wigner 
resonance shape for the top quark, we obtain an infrared finite correction to the 
mass distribution. The major effect of the QCD radiation is to degrade the peak, 
with practically no shift in the position of the maximum. For energies not too far 
above the tt threshold, most of the gluon radiation occurs during the decay of the 
quarks; however, at higher energies the radiation off the tops during the production 
phase becomes more important. 
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Appendices 

A e+e~ tig production amplitudes 

The real radiative corrections to ti production and decay can be given by helicity 
amplitudes, with only minor complications due to the three-body final state. We 
can describe the ttg production event configuration in the lab frame in terms of five 
variables. Two of these are the energy fractions Xi = 2Ei/ ^/s of the top and the 
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gluon, which are in turn determined by the variables of equation (|T8|). These fix 
all of the lab-frame energies and angles within the tig-p\a.ne. Two more variables 
are just the polar angle 6 and azimuthal angle of the top quark with respect to 
the electron beam axis. The final variable that we need is the angle (pg between 
the e'^'e^t-plane and the ttg-p\ane, rotated around the top momentum axis. Note 
that the rotation by (pg around the top quark momentum axis also rotates its decay 
products. This completely determines the event kinematics. 

For longitudinally polarized electrons, the intermediate photon-Z state will be 
an eigenstate of spin along the beam axis. However, it is more convenient to work 
in a basis where the vector boson is a spin eigenstate along the top momentum 
direction. Labeling these eigenstates by 7a, we can expand the matrix elements in 
terms of amplitudes in the new basis, which are now independent of the variables 
(j), 9, and (pg-. 

'T^lM{L- 7i ^ ttg) + T^rM{R; 7l ^ ttg)] ^(1 + cos^)e 
y^LM{L- 7« ^ ttg) + r^RM{R; 7ij ^ tig)] -^(1 - cos 0^^'^^ 
'j^tiML; iz ttg) + 7z ^ ttg) 

MieRCL ttg) = e' 

- y^LM{L; 7L ^ ttg) + J^^rM^R; 7l ^ tig) 

- 'j'^LMiL- 7^ ^ tig) + J'^RMiR; 7fl ^ tig) —(^^ 
+ f^ilA1(L; 7z tig) + ^ilA^(i?; 7z ^ tig) 



+ 
+ 



sin^ 



-i(j>g 



(2i 



71' 



cos( 



-^(1 + cos, 



sin^ 
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We have also separated the pieces arising from the left-handed and right-handed 
currents. The form factors = T\y + T\j^ and T\i^ = T\y — T\j^ are obtained 
from equation (H) evaluated at tree level. 

The matrix elements in equation (^) with left-handed currents are: 

. ^tg 6tg dtt f n , /i \\ 

sm ^ COS ^ COS— [xtPt + (1 -Xt)) 

sm — sm — [1 — Xt) 
2 2 

sm —{xtPt COS — sm Y + (1 - ^t) sm ) 

sm —{xtPt cos — cos + (1 - Xt) cos ) 





tLhgL) = 




M{L;-fn - 


tLhgL) = 




M{L;-fz - 


tLhgL) = 


V2 


M{L;^L- 


tLhgn) = 




Mm-fR - 


tLhgn) = 





M{L;-fz - 


tLhgn) = 


A 

A+- 
V2 


M{L--fL- 


tnhgL) = 


-A__ 


Mm^R - 


tnhgL) = 




M{L;-fz - 


tnhgL) = 


V2 


M{L-^L- 


^ tRtLgn) = 





M{L;-fR - 


tRhgn) = 




M{L;-fz - 


tRhgR) = 


A__ 
V2 


M{L;^L- 


tLtRgL) = 




M{L;-fR - 


tLhgL) = 




M{L;-fz - 


tLhgL) = 


72 



ett 


- Otg 




2 


ett 


+ ^tg 




2 


Ott 


+ Otg 



dtg f n ■ ^tg ■ ^tt . ,^ _ n 

COS Y-l^^tPi sm Y"Sin Y + (1 - Xt) COS ^ 

2 dtg Ott X 
COS —COS [1 — Xt) 

2 2 ^ ^ 

. Otg 6tg . Ott/- n /I \\ 

sm — COS — sm—[xtPt - (1 - Xt)) 

Otn I ^ . Ota Ott /-, \ . Ott — Ota-. 

COS —[XtPt sm Y cos Y + (1 - Xt) sm ^ ^ 



Otn , ^ . Ota . Ott / — \ Ott ~l~ Otq , 

COS —{xtPt sm Y sm Y + (1 - a^t) cos , 

■ ^ig t n Otg Ott , _ s Ott + Otg 

sm — [XtPt cos — cos -— + (1 — Xt) cos 



2 ^ ' 2 2 ' ' 2 ' 

(29) 



%g Otg . Ott 



2 2 2 ' ' ' 

. 2 ^tg Ott /-I X 

sm —COS — [\ — Xt) 
2 2 ^ ' 

Otq I n Ota Ott /-, \ Ott — Otq^ 

sm — {XtPt COS — COS 1- (1 — Xt) cos 
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M[L]'^L ^ thtngB) = A++sin — (xt/^j cos — sm Y + (1 - xt) sm ) 

Mm-fR-^tLtRgn) = 

M{L]-iz -^tLtRgR) = — -j=-cos—{xtfJtsm — cosY-{^-Xt)sm ) 

M{L;'jL tntRgL) = cos^ ^ sin y (1 - x*) 

M{L;^R^ tRtRgi) = A_+ silly cos y cos - (1 -Xt)) 

M{L]-iz ^tRtRgL) = — -y cos-^(xt/:/tSin^smy - (1 -xjcos — ^) 

M{L;-fL-^tRtRgR) = 

M{L;-fR-^ tRtRgR) = -v4_+ cos y(xt/:/t silly cos Y - (1 -xjsm ) 

M{L;'jz ^ tRtRgR) = — y silly (xt/^/t COS y silly + (1 -xtjsm ) 



where 



Xg 



with Tr (T°T ) = 6 /2. The remaining matrix elements can be obtained from 

M{L, R- 7a ^ tLtL(?,) = -(-l)^A^(i?, L; 7_a ^ tfitR^?-.) 

7a ^ tLtfl^.) = (-l)^A^(i?, L- 7_A ^ tfltL^?-.) • (31) 



In terms of the variables in equation (p^ ) the energy fractions are 



2 ■ 2^V 1-x 

Xf 2 Xq- 



= l-7f + ^.(^-i) (32) 



"9 



Here = 1 — 4m^ /s is the tree level velocity of the top quarks, while the velocities 
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of the t and t in the presence of the radiated gluon are 

4m? 



= 1- 



(51 = (33) 



The lab-frame angles are obtained from 



1 r - - '^^t^ 
COS ^tt = o - o - Xt- Xt + XtXt H 

COS 6tg = — ^ — [xt - xt- Xg + XtXg] . (34) 



B t ^ bW^g decay amplitudes 



The helicity amplitudes for top decay with a radiated gluon can be calculated in 
an analogous manner to the production calculation in appendix A. We describe the 
decay configuration in the top rest frame in terms of five variables. Two of these 
are the energy fractions Xi = 2Ei/mt of the W'^ and the gluon. These energies 
are determined by the variables of equation (|21|), and they fix all the energies and 
angles within the hW^g decay plane. Two more variables are the polar angle Xt 
and azimuthal angle ipt of the with respect to the top momentum boost axis. 
The final variable is the angle (pg between the plane given by the top boost axis 
and the momentum and the 6iy"'"(y'-plane, rotated around the momentum. 
This rotation by (pg also rotates the decay products. 

We can make explicit the dependence of the matrix elements on the variables 
Xt, 'ipt^ and (pg if we expand the top helicity eigenstates th onto a basis of spin 
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eigenstates along the momentum direction. Labeling these states as t'^ we 



obtain the relations: 




MitL- 


bW+g) 






M{t'R 


bW+g) 


MitR- 


bW+g) 






M(4 


^ bW^g) 



, ^e^'f'o/^ + M{t'L^ bW+g) cos ^e-^<^«/' 



(35) 



. ^e^0./2 _ _yv4(i^ ^ bW+g) sin ^g-^^^/^ 
The helicity amplitudes in this basis are 



Mit'R 
M{t'R 



M{t', 



bLWngL) 

bLWngL) 
bLWigi) 
bLWLgL) 



M{t'R^bLWzgL) 
M{t'^ ^ btWzgL) 



bhWRgR) 
biWigR) 
bLWign) 



M(4 ^ bLWzgR) 



2 / 9wg , Owb ^Wg 
1= Xn cos + Xh COS COS — - 

' 2 2 












( . Qwg , . Qwb Owg + ^VF6\ 

— sm — - + Xft sm cos — 



\f0^g 

Xw{^+(^w)( . ^Wg , . 

— ^— — y-Xg sm + Xfo sm — 



W^2CXg 



9wb ^Wg 

COS — - 



Owg + Owb 



2 

— / dwg + (^Wb\ 
; [Xg COS — ^ + Xfe COS — - COS ^— ) 







6*1^6 / a;^ + 6wb . ^ \ 

COS — Smt^vKg ) 



/ — ^Wb^-. a \ 

-Va;6C0s^-(l - co^dwg) 
-i/xfosm ^-(1 + COS&'iyg) 



(36) 



-^/i^sin 2 



2\/2i 



0Wb /„ / a;;, 9wg + ^W'b , • /I \ 

— y \l~ ''^^ — 9 — + ^'"^ i 

1 + Msm— (- 



^H/fe / o / 6wg 
COS — — 



(^Wb 



+ sin 6wg^ 
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+ (1 - Pw) COS ^^(1 + COsOwg] 



2^2- 



w 



{I - fJw) COS ^[2^1 — COS l-smOwg 

0Wb , 



+ (1 + Pw) sin ~ cosOwg 

where we have dropped a factor of —iT^-g^g / V2. In terms of the variables of (|2 
the energy fractions are 



Xg = y{l-w' 

1 I 2 ^g{^ — W — Xg) 

Xw = I + W — Z — — (37) 

1 — Xg 



and we have also introduced the variable ( = 2pb ■ Pg/mf = 1 + w"^ — xw The 
velocity of the is given by 

Pi = 1 - ^ . (38) 



The angles in the top rest frame are obtained from 

cos 9wb = n i^q - Xw - Xb + xwXb + 2nP'] 

XwPwXb 

I 

cos 9wg = n [Xb — Xw — Xg + XwXg + 2w'^] . (39) 

XwPwXg 

The amplitudes for i decay in its rest frame can be obtained from these by 
simply using 

M{th-^b,W^g^) = Mii^h^b-.Wl^g-a) , (40) 

while replacing all of the energies and polar angles of t decay with the corresponding 
variables of i decay and replacing the azimuthal angles by tpt —'ipt and (pg ~4>g- 
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/ic«t (in GeV) 
5 10 20 oo 


production 
only 




-2.8 


-0.7 


+1.2 


+3.3 


+3.2 




+4.2 


+3.6 


+3.0 


+2.5 


+2.9 



Table 1: Percentage 0{as) corrections to the top quark forward-backward asym- 
metry for rrit = 175 GeV and y/s = 400 GeV with polarized electrons. The first 
four columns are using the reconstruction algorithm (|2^) , while the last column 
gives the corrections from production only, assuming exact event reconstruction. 
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2: Polar angle dependence of from decay of right-handed top quark. 
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3: Polar angle dependence of from decay of in a tt event. 
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Fig. 4: Fraction of events containing a production gluon as a function of ^/s. The 
curves are, from top to bottom, for /icut — 5, 10, 20, 40, 80, and 160 GeV. 
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Fig. 5: Top mass reconstruction distributions for ^/s = 400 GeV (a) in the zero- 
width hmit and (b) with an initial Breit-Wigner resonance distribution. The his- 
tograms are for ficut ~ 5 GeV (dots), 10 GeV (dotdash), 20 GeV (dashes), and oo 
(sohd). The smooth curve in (b) is the original Breit-Wigner distribution. 
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Fig. 6: Top mass reconstruction distributions for ^/s = 1 TeV with an initial Breit- 
Wigner resonance distribution. The histograms are for j^cut = 5 GeV (dots), 20 
GeV (dotdash), 80 GeV (dashes), and oo (sohd). The smooth curve is the original 
Breit-Wigner distribution. 
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Fig. 7: The 0{as) corrections to the top quark polar angle distributions for y/s = 
400 GeV with (a) left-polarized electrons and (b) right-polarized electrons. The 
histograms are for fi^ut = 5 GeV (dots), 10 GeV (dotdash), 20 GeV (dashes), and 
oo (solid), while the points plotted with the symbol x are the pure production 
corrections. 
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Fig. 8: Distribution of observed cosxt for events with true cosxt between -0.1 
and 0.0. The histograms are for Hcut = 5 GeV (dots), 10 GeV (dotdash), 20 GeV 
(dashes), and oo (sohd). 
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Fig. 9: Top mass reconstruction distributions without energy smearing of the final- 
state partons for ^/s = 400 GeV (a) in the all-jet mode and (b) in the lepton-|-jets 
mode. In both plots the dotted histogram is at tree-level, the solid histogram is at 
0{as), and the smooth curve is the original Breit-Wigner distribution. 
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Fig. 10: Top mass reconstruction distributions with energy smearing of the final- 
state partons for y/s = 400 GeV (a) in the all-jet mode and (b) in the lepton-|-jets 
mode. In both plots the dotted histogram is at tree-level and the solid histogram 
is at 0{as). 
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